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CONDITIONS OF REACTIVE-DIFFUSIVE TYPE 
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Abstract. This paper deals with the heat equation posed in a bounded reg- 
ular domain f! of (TV > 2) coupled with a dynamical boundary condition 
of reactive-diffusive type. In particular we study the problem 

{ut — Am = in (0, oo) X SI, 

ut = kun + /Apn on (0, oo) X F, 

m(0, x) = uo{x) on r, 

where u = u{t,x), t > 0, x S O, F = dQ, A = A^; denotes the Laplacian 
operator with respect to the space variable, while Ap denotes the Laplace— 
Beltrami operator on F, v is the outward normal to n, and k and I are given 
real constants, I > 0. Well— posedness is proved for data mq G H^{Q) such that 
'^Olr G ^^(r). Wo also study higher regularity of the solution. 



1. Introduction and main results 

We deal with the evolution problem consisting in the standard heat equation posed 
in a bounded domain, supplied with a dynamical (or Wentzell) boundary condition. 
The precise problem is 

{Ut — Au = in (0, oo) x Q, 

Ut = kuu + /Apw on (0, oo) x F, 

u(0,2:) = uo{x) on 57. 

Here u = u{t,x), t > 0, x G fl, where 51 is a C°° regular bounded domain of 
(iV > 2) and F = 957. The first equation states the law of standard diffusion or 
heat conduction in 57, and A = denotes the Laplacian operator with respect 
to the space variable. In the boundary equation the value of u is assumed to 
be the trace of the function u defined for x G 57, Ar denotes the Laplace-Beltrami 
operator on F, i/ is the outward normal to 57, and fc S M and I > are given 
constants; the term ku^ represents the interaction domain-boundary, while lA^u 
stands for a boundary diffusion. 

A number of authors have studied parabolic problems with dynamical boundary 
conditions like Note that we can replace ut by Au in this boundary condition 
which leads to the form known as generalized Wentzell boundary condition. The 
problem has been mostly studied the case when there is no Laplacian term on the 
boundary condition, i.e., when 1 = 0. In particular, when fc < Problem ([1]) is 
well-posed. See [l], [12], [15], [l6], [17], [H], [IS], [23] [24] in the case fc < which 
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represents a dissipative interaction; the non-interactive case /c = is rather trivial. 
However, when k > we are in the presence of a reactive interaction and Problem 
([T]) is ill-posed, as shown in the recent papers |3] and [35]. See also [5] and [21] for 

the related case k = k{x). 

The question we address in this paper is the following one: is the situation im- 
proved by adding to the dynamical boundary condition a Laplace-Beltrami correc- 
tion term with I > 0? The interest of such a correction both for the modeling of 
parabolic and hyperbolic problems has been recently pointed out in [22] . In par- 
ticular ([ll describes (see [221 P- 465]) a heat conduction process in with a heat 
source on the boundary which can depend on the heat flux around the boundary 
and on the heat flux across it. The case of dissipative interaction, k < 0, has been 
studied in [3], [TO], [20] (see also [7], [5] and 123 )■ It turns out from the quoted pa- 
pers that Problem ([1]) is well-posed in the framework of LP(yt) x i^'(r), 1 < p < oo. 
This is to be expected since both terms in the right-hand side of the boundary 
condition have the "favorable sign" . The aim of this paper is to solve the system 
in the reactive case A: > 0, that is in the usually ill-posed case. The estimates of 
the quoted papers did not allow to cover this case. 

A first step in this study has been performed by the authors of the present paper 
in [33], where we consider the Laplace equation instead of the heat equation as 
domain equation. The modified problem admits a simple functional framework; 
the paper helped the authors understand the dynamical boundary condition (|T])2 
and allowed us to formulate the conjecture that turns out to be correct, but the 
arguments used there do not work for the heat equation. Indeed, a new estimate 
is needed to deal with Problem ([1]), which cannot be obtained in the framework of 
L^in) X L^{T). 

We want to show that Problem (|T]) is well-posed in an appropriate setting. We 
propose to work in the space 

(2) H = {{u,v)eH\n)xH\T):u^r = v}, 

where U|r denotes the trace of u on F, with the natural topology inherited by 
H^{n) X H^{T). Here and the sequel, we denote for any s e M, H'in) and i7"(r) 
the Sobolev spaces of complex-valued distributions respectively on D, and F (see 
[26] or [30]). For the sake of simplicity we shall identify, when useful, H with 
is isomorphic counterpart {u E H^{n) : u\r <= H^(r)} through the identification 
(u, U|r) I— >■ u, so we shall write, without further mention, u Cz H for hmctions 
defined on Q. 

Our main result is the following 

Theorem 1. For any uq €z H Problem ([T]) has a unique solution u = u(uq) such 
that 

u e C([0, c5o); H^{n)) n Ci((0, c5o); H^{Vi)) n C((0, c3o); ff^(f^)) , 
u\r e C([0, oo); i?i(F)) n C\{Q, oo); H\T)) n C((0, oo); H^(T)) . 
Moreover, 

(4) l|Vti(<)|li.(f,) + Wdru^rimhiT) + \\u\rml^r) 

< e'^"*l|Vuo|li2(f2) + \\druo\r\\h^r) + ho|rlli2(r) 
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for all t > 0, where Xo > is a constant depending on D,. Finally, the family of 
maps {uo I— > u(ua)(t),t > 0} extends to an analytic quasi- contractive semigroup in 
H , and consequently 

(5) M e C°°((0,oo) X H). 

The solutions are in principle complex-valued but it is clear that for real-valued 
data the solution is likewise real-valued. As usual, more regular solutions are ob- 
tained for more regular initial data satisfying usual compatibility conditions. This 
is the content of the following regularity result. 

Theorem 2. If uq G H'^"'+^{n) and uo|r G i?^"+^(r) for some n S N, and 

(6) (A*Mo)|r = fc(A*"'uo). + ;Ar((A'-iuo)|r), for all i = 1, . . . ,n, 
then 

u e C([0, oo); n Ci([0, oo); H'^'''-\n)) n . . . n C"([0, oo); H^{n)), 

u\T € c([o, oo); i?2»+\r)) n c\[Q, oo); iJ^'^-^r)) n . . . n c"([o, oo)- h\t)). 

Finally, if uq G C°°(n) and ^ hold for all i £ N, then 
(8) uG C°°([0,oo) xH). 

The proofs of Theorems [T] and [2] rely on the study of the resolvent problem with 
eigenvalue-dependent boundary condition, that is 

I — Au + \u = h in J7, 

1 —kui, — ZAru + \u = h on F,. 

where A G C and h £ H. Such type of problems has been studied by some authors, 
starting from the classical papers (see [13], [M]) to more recent ones (see [4] and 
the bibliography therein). Our result concerning Problem ^ is Theorem [3] below. 
Finally, we study the limit behavior of the solution u when I — > 0+ (vanishing 
boundary dissipation). See Theorem [HI below. 

The paper is organized as follows. In Section [5] we recall some well-known facts 
and we state some preliminaries. In Section [3] we analyze the elliptic problem 
while in Section 2] we apply the results obtained to Problem ([1]). In Section [5] 
we analyze the limit behavior when ^ ^ 0"*", while the final section contains some 
comments on future developments. 



2. Preliminaries and functional setting 

Notation. We denote by || • ||p, 1 < p < oo, the norm in L'''{il) and, also the norm 
in L''(ri;R^) since no confusion is expected. We denote by || • ||p_r the norm in 
LP{T) and also, when p = 2, the norm for square intcgrable 1-forms on F. 

Laplace— Beltrami operator. We recall here, for the reader's convenience, some 
well-known facts on the Laplace-Beltrami operator Ap. Wc refer to [5S] or [30] 
for more details and proofs. We start by fixing some notation. Clearly, F is a 
Riemannian manifold endowed with the natural metric inherited from , given in 
local coordinates by {gij)ij=i^...,N-i- We denote by dV the natural volume element 
on F, given in local coordinates by y/g dyi . . . dyN-i, where g = det{gij). We denote 
by Vr the Riemannian gradient and by dr the total differential on F. We use the 
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notation (•, •) for the Riemannian inner product of vectors while (-j-) is used for the 
natural scalar product on 1-forms on F associated to the metric. Then, it is clear 
that {dru\drv) — (Vru, Wrv) for u,v G C^(r), so the use of vectors or forms in the 
sequel is optional. 

The Laplace-Beltrami operator Ar can be at first defined on C°° (F) by the formula 
(10) - J {Aru)vdV ^ J {dru\drv)dV 

for any u,v E C°°{T), and it is given in local coordinates by 

where (g*^) = {gij)~^ as usual. Clearly, by pT|) . Ar can be considered as a bounded 
operator from H'^^^{T) to H^{T), for any s e M. Consequently, Formula (fTO)) 
extends by density to u,v G if^(F), where the integral in the left-hand side has to 
be interpreted in the distributional sense, as Aru G H~^{r). 

Remark. In the sequel, the notation dV will be dropped from the boundary 
integrals; we hope that the reader will be able to put in the appropriate integration 
elements in all formulas. 

Since Apl = the operator is not injective, but by (|T0| we have 




so that the operator L := — Ar + 1 is a topological and algebraic isomorphism 
between H^{r) and if^^(F). Moreover, by elliptic regularity (see [301 P- 309]), 
i-i : F'=-i(r) ^ H''+\T), fc = 0,1,2,..., is bounded, so L : iJ'=+i(F) ^ H'^-^T) 
is an isomorphism. By interpolation, : H'^{r) H'''^^{T) for all s G M, s > — 1, 
giving the inverse of L : iJ^+^(r) i/*(F). By duality, this fact holds for all real 
s. 



Dirichlet— to— Neumann operator. We will also need some well-known facts 
about this operator that will be used at some technical points. We refer to [55] 
for details and proofs. For any u E i/*(F), s € M, the non-homogeneous Dirichlet 
problem 

(13) 

I w = u on F, 

has a unique solution v € H'^~^^^'^(fl), here denoted hy v = Dm. Moreover D is a 
bounded operator from _ff''(F) to _ff''+^/^(il) for all real s, and has a normal de- 
rivative Vi, £ H^~^{T). The operator u i— > v^, known as the Dirichlet-to-Neumann 
operator, is bounded from H'^{r) to 7?^~^(F), and it will be denoted in the sequel 
by A. For all u,v £ C°°{T), integrating by parts twice we have 

(14) Auv^ / V(Du)V(Di;) = / uAv 

Jr Jn Jr 



which, by density, holds for all u,v £ iJ^(F). 
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Functional setting. In the sequel we equip H^{T) with the equivalent norm in 
(fT2l). so we denote 



(15) 



{u,v)Hi{r)^ / UV+ / {dru\drv), = (u, -u)//i(r) 



for all u,v G if^(r). Moreover, since — Ar + 1 : H^{r) i^(r) is an isomorphism 
we can equip H'^{T) with the equivalent norm 



(16) 



(r) 



uv - 



for all u, V E H^{T). Moreover, we denote as usual 



(17) 



(Q) 



The SPACE H. We now introduce, as anticipated in the introduction, the space H 
given in ([2]) , which by the Trace Theorem is a closed subset oi (fi) xH^{T), hence 
a Hilbert space with respect to the scalar product inherited from H^{n) x H^{r). 
For the sake of simplicity, we shall drop the notation u\y, when clear, so we shall 
write ||m||2Tj /r ^^'^ so on, for elements of H, through the already mentioned 
identification (m, U|r) u. We equip H with an equivalent norm which simplifies 
our calculations. This is the content of the following 

Lemma 1. We set, for any u,v £ H, 



(18) {u,v)h 



VuVv ■ 



{dru\drv) 



(u, u)h- 



The 



\h is equivalent in H to the standard norm inherited by H^(Q,) x H^{T). 



Proof. We just have to show that if we drop 1 1 • 1 1 2 in the standard norm of ^ ($7) x 
H^(T) we get an equivalent norm. This follows by a Poincare-type inequality which 
says (see [Ml Theorem 4.4.6] in the real valued case, the extension to the complex- 
valued one being trivial) that 

<Ci||Vu||2 for all w e iJ^(17), 

where Ci = Ci{N, n) > 0, 2* is the Sobolev critical exponent, i.e. 2* = 2N/{N-2) 
when iV > 3, 1 < 2* < 00 when N ~ 2. Consequently, since il is bounded and F is 
compact, we get 



u — 









(19) 



\H\2< 



<Ci\\Vu\\2 + XN{^) / |"| 



<C2(||Vlij|2 + ||u||2,r) 

where A^r denotes the usual Lebesgue measure in and C2 = C2{N, U) > 0. This 
estimate completes the proof. □ 



The space V. We need a further space 

(20) V = {{u,v) £ H^{n) X H'^iV) -.u^r =v} 

which is naturally embedded in and it is a Hilbert space with respect to the 
scalar product and norm inherited from H^{n) x H^{T). As before we equip it with 
a suitable scalar product which induces a norm equivalent to that one. 
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Lemma 2. // we set, for any u, v e V , 
(21) {u,v)v 



/ AuAv + 


J ApuArv + 




\\u\\y {U,U)V 


Jn 









then II • ||y is equivalent in V to the standard norm inherited by H^{il) x H^{T). 

Proof. It simply follows by elliptic regularity estimates. Indeed, for any u £ H^{il) 
we have (see [H p. 202]) 

\\u\\H^n) < C3 {\\Au\\2 + ||u|rllH3/2(r)) 

and consequently, since H^{T) is continuously embedded in H^/^{T), for any u <£V 
we get 

(22) WAhhu) < \\Au\\2 + \\u\\H2(r) 

which by ([T5| completes the proof. □ 

3. Elliptic theory 

This section is devoted to study the solvability of the coupled elliptic system ^ 
when I > 0, k e R, X e C and h e H . 

Definition. By a solution of Problem ^ we mean a function u £ V such that ©I 
holds true in L'^{ft), while ^2 holds true in L^(r). 

Space V was just introduced in (|20p . Before stating the main result of this section 
we introduce, for any s > 1, the further space 

(23) = {{u, v) e H'in) X iJ"(r) : w|r = v}. 

Clearly, being closed in the product space i?*(i7) x i?*(r), H'^ is a Hilbert space 
equipped with the norm inherited norm, which we denote by || • ||^fs. Moreover, 
it is naturally embedded in H and = H, ~ V (more precisely, || • ||^i and 
II • \\h are merely equivalent, like || • ||/f2 and || • \\v)- 

Our result concerning Q is the following 

Theorem 3. There is a positive constant Aq, depending on /, fc, fl, N , such that for 
A G C, ReX > Ao and any h € H Problem ^ has a unique solution u Cz V , which 
also belongs to . Moreover, if h £ iJ" for some s > 1, then u G H"^^ . 

Finally, there is C4 ~ C^ll, fc, f2, s. A) > such that 

(24) ||M||a.+2 < C4||M||ff= forallheH". 

In order to solve elliptic problems via the variational method it is useful to intro- 
duce a sesquilinear form, which leads to weak solutions. The most natural way to 
perform this procedure for Problem ([9|) would be to multiply (at least formally) the 
equation — Ait + Au = /i by a test function (jj and integrate over to get 



Aucj) + A / u(j) = / h(\). 
n Jn Jn 

Integrating by parts, when u is regular enough, 

VuVu — / u,^<f) + X u(j) = / h<j>. 
n Jr Jn Jn 
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Then, using the boundary equation in ([9]) we get (when fc 7^ 0) 

VuVv + Y [ hlj)- y [ u(i)+ y [ Arucj) + X [ ulj) = I hcj). 



ft k Jy k Jy k JY Jn 

Finally, by PH)) we arrive to 

(25) / VuV0-| / (drw|<ir0)-T ( u4> + X ( u4> = ( hlj) + ( hlj). 
Jn k Jy k Jy Jn k Jy Jn 

Now, it is easy to check that the sesquilinear form in the left-hand side of ([^5)) is 
indefinite in the case fc > 0, so this procedure does not produce useful estimates. 
Thus, one has to look for a positive definite sesquilinear form, at least for ReA large 
enough. This is exactly the content of the following two lemmas. The first one 
introduces the sesquilinear form which turns out to be appropriate. 

Lemma 3. Let h e H . Then u Cz V solves Problem ^ if and only if 

(26) a\{u,v) = {h,v)H for all v &V, 
where the sesquilinear form ax on V is defined by the formula 



(27) 



ax{u,v)^ / AuAv + l / AyuAyv + X / VuVv + {X + l) / (drwldrw) 
Jo Jr Jn Jr 

—I / Aruiv + fc / u,^ Ayv — k / u^v^—k Ui,v + Xl uv. 
Jr Jr Jr Jr Jr 



Moreover in this case u G H^(fl) and 7i|p E H^(T). 
Proof. It is divided into several steps. 

(i) Claim. If u e V is a solution of ^ , then u e {il) and u\y {^) ■ To 
recognize that our claim is true we use elliptic regularity both on and F as follows. 
Since u S H^{n) we have u^, € H^/^{T) by the Trace Theorem. So, being H^y G 
H^{r) and e i/^/^(F), from ^2 it follows that -Ayu + u e H^/'^{T), so that 
using the the isomorphism property of —Ay + 1, we conclude that U|r G 77^/^ (F). 
Consequently, using elliptic regularity for nonhomogeneous Dirichlet problems f |26l 
p. 203]) we obtain by that u e H^{n). 

From this, and using the Trace Theorem again, we get <E H^^^{il.). Using ©2 
again we then get —Ayu + u E H^(r), so as before U|r G H^(r), completing the 
proof of our first claim. 

(ii) Claim. If u E V is a solution of then formula ([^5]) holds. By the first 
claim we have Au € H^{^1). Moreover, by ©I we get 

(28) (AM)|r = AM|r-/i|r ei7^(F). 
Consequently, we get that Au G H, so from ©i we have 

(29) {-Au,v)h + X{u,v)h = {h,v)H for all i> € i7. 
Formula (|29|) can be written more explicitly, using (|T8)) . as 

(30) [ V(-Au)Vd+ [ {dY(-Au)\dYv) - [ Auv 
Jn Jr Jr 

+ X VuVu + A / {dYu\dYv) + X / uv^{h,v)H. 
Jn Jr Jr 



8 



JUAN LUIS VAZQUEZ AND ENZO VITILLARO 



Now, using ([28| we can write ^2 in the form 

(31) {Au)\r ^ ku^ + lAru\r- 

Plugging ([3T|) into (|30l) we get 



(32) / V(-Au)Vw-fcy {dru^\dTv)-l J {drATu)\drv) - k J u^v 

-I / Atuv + X / VmVij + A {dru\drv) + X / uv ^ {h,v)H 



for all V Cz H. Now we restrict to test functions v e V, we integrate by parts the 
first integral in (|32l) and we use (fTO|) in the first one to get 



(33) / AuAf — / Auvi, — k / (dru,y|drw) + ^ / ApuArv ^ k / u^w 
Jn Jr Jr Jr Jr 

-I / Aruw + A / VuVw + A / (drwldrw) + A / uv^{h,v)H- 
Jr Jn Jr Jr 

Plugging ([5T|) once again in the second integral in the left-hand side of and 
([TU)) in the third and sixth ones we finally get (pS)) . 

(iii) To complete the proof, we now suppose that (|26)) holds for some u £ V. We 
have to prove that u solves (|9]). Integrating by parts in the third integral in ((27)) 
and in the first one in (1181) we can then write (l26l) as 



/ AuAv + l / ArwAru-A / uAv + X / utJ;7+(A + {dru\drv) 
Jn Jr Jn Jr Jr 

(34) —I / Aruw^+A; / u^Ayv — k i UyV^ — k i u^v + X / uv 



n 



hAv + I hv^+ I {drh\drv) + / hv 

for all w G V. Using ([TU)) we can write (|M)) as 

AuAv + I / AruArw — A / uAu + A / uTv — X j uArv 
n Jr Jn 



n 



(35) — ^ / Ayuv — I I ApfiWi/ + fc / u^Ayv — k j u^v^ — k J Ui,v + X J uv 

hAv + / /iTJT- / /iAr?J+ / /iiJ, 
Jr Jr Jr 

that is, by grouping the terms with respect to the test function. 



(36) / {Au-Xu + h)Av+ / {-I Aru - ku^ + Xu - h) {-Arv + v^ + v) ^ Q. 



The form of ()36)) suggests now how to proceed. Indeed if we restrict to test 
functions v G C^(il), at least to get we get that A(— Ait + Xu — h) = in 
distributional sense, which is not Then it is more useful to start by proving 

^2- With this aim, we restrict (|55)) to test functions EDw, where v € H'^iT), to get 

(37) j {-IAtu - ku^ + Xu- h) {~Arv + Av + v) = 0, 

where A denotes the Dirichlet-to-Neumann operator already introduced. 
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We now claim that by ([ST]) it follows that 

(38) {-lATU~ku^ + Xu-h)'$=0 ioi al\(j>e L\r), 



from which clearly one has that ^2 holds in i^(r). To prove our claim it is enough 
to recognize that, given an arbitrary (j) G L^(r), the problem 

(39) - Arv + Av + v^(j) 

has a solution w G H^(T), which turns out to be unique. Hence, our claim is 
nothing but a refinement, in this particular case, of a previous result of the authors 
[551 Lemma 1] which says that given I > and /c G K there is A > such that for 
A > A the problem 

(40) - I Arv - kAv + Av = (j), 

has unique solution v G H^{T) for any ip G i^(r). In particular, our claim is proved 
if we prove that, when fc < 0, then we can take A = 1. To prove this fact we argue 
as in the quoted paper, writing (|40p in the more explicit form 



(41) - J kAv^J + I j {dvv\dTtp) + K j v^j ^ j (t)ij for all ^ € (F), 
and then we apply Lax-Milgram theorem (see fill p. 376]) to the scsquilincar form 
a(w,'0) = - / kAv4) + l / (dri'Mr'0) + A / v4), v.ip € H^i^), 



which is trivially Hermitian (by (jl4p ) and continuous. To recognize that it is also 
coercive for A > 1 we simplify the argument of [33] . Indeed, since fc < we have 
by da 

a{v,v) = -fc||V(Dij)||2 + [||dri'||2,r + A||i'||2,r > min{f, A}||i;||^i(p), 

so that the form is coercive whenever A > 0. Then, from Lax-Milgram theorem we 
get the existence of a solution v G H^{T) of (|40p . By the isomorphism property of 
— Ar + 1 is then follows that v G H^{T), completing the proof of our claim. 

Now, to prove we use ^2 in (|36)) to get 

(42) / (Au - Xu + h) Av for aU v eV, 

Jn 

which clearly implies since for any ip G L'^{n) there are v E V such that 
Av = ip, for example by taking the unique solution v G H'^{VL) n Hq{VI) with 
homogeneous Dirichlet boundary conditions. The proof is then complete. □ 

The following key estimate shows that the scsquilincar form (|27p is appropriate. 

Lemma 4. There are positive constants Aq and C5, depending on l,k,n,N, such 
that for all X G C, ReX > Xq we have 

Reax{u,u) > C^WuWy forallueV. 
Proof. By ((27|) we have 



axiu, u) =11 Aii||^ + /|| Arfilllr + A||Vu||^ + (A + l)\\dTu\\lj^ + X\\u\\l^-p 
^l / Auu^ + k / UjyAru — k\\ui^\\'2 -p — k / u,/u 



r 
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and then 

ReaA(M,u) ==||Au||2 + ;||Aru||2,r + R-eA||Vu||2 + (ReA + /)||dru||2,r 

+ReX\\u\\l Y + {k - I) / Re[AruUi,] - /e||w^||2.p - fc / Rc[M„It] 

(43) ^ ^ 
>|1 Au||^ + ;|| Arulla.r + R-eAII Vu||2 + (RcA + Ollrfritlllr 

+ReA||u||^_r - (1^1 + J \^ru\\u^\ - \k\\\u^\\lr ~ \k\ J\u^\\u\. 
By Young inequality we estimate 

(44) |fc|^KII«l< Yll«.ll^,r + Ylkll2,n 

and, given any e > to be fixed later, by weighted Young inequality 

(45) i\k\+l) |^|Ar.||..|<M±ll£||AH|^^r + ^h.||,r. 

Plugging dH]) and (05]) into (gS]), we get 

i\k\+l)e] 



ReaA(u, u) >|| Aw 



2 



2 



||Arw||ir+ReA||Vuii2 



(46) + (rcX \\u\\l^ + (RcA + mdruWl, 



+ 1 3 



1^1 h.llir- 



2e 2 

Then, by choosing e = l/{\k\ + /), we get 

^^^^ ReaA(u,w)>||Au|!2 + l||A^^t||l^ + RcA||V^i||2+ (^RcA-i|) ||u|||r 
+ (ReA + /)||dru||2_r-C6||M,|l^,r, 

where Cg = Ce{k,l) ~ (■^^^^2r^ — ^ ll^'^l)- To estimate the last term in the right- 
hand side of (|T7|) . we note that by the embedding H'^/^{Q) ^ H^^'^{ri) and by the 
Trace Theorem there is C7 = C7(f2) > such that 

||w^||2,r < C'7||w||^7/4(f2) for all u e H'^{n). 

Consequently, by interpolation inequality (see |26]). 

Using weighted Young inequality we then get, for any 5 > Q (to be fixed below), 
lk.|l^,r < j^Mm(n) + ^^Ml^^n) for all u e H\n). 



By applying ([T9| and ([22]) in the last formula, we get 
5 



(48) h.||^,r < ^(IIVt^II^ + \\u\\l^) + C8<5 (||Ati||2 + \\^ru\\l^ + l|«||^,r) 
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for all u eV, where Cs = Cs{^) > 0. Pluggmg 

I 



into (l47l) we derive 



(49) 



Rea^{u,u) >{l ~ CeC86)\\Au\\l + -(2 
+ {ReX-C6Cs/5)\\Vu\\l 



{ReX + l)\\dru\\ 



2 

2,r 



ReA - 



CeCs (6 + 1/6) \\u\\Ij 



Choosing 6 = 6o = min {2, 1} / {AC^Cs) we rewrite as 

ReaA(u,M) >i||Au||2 + ^||Aru||^ r + (ReA - CeCs/6o)\\Vu\\l. 

(50) 2 4 

+ (ReA + l)\\dru\\lj. + \rcX - ^ - CgCg (<5o + l/(5o) 

Now, by setting 

(51) Ao = max {CeCs/6o, \k\ + 2CeCs (^o + V-^o)}, 
we clearly have, when ReA > Aq, that 

|fc| 



2,r- 



RcA - CgCg/^o > and RcA 
so by ([50| we finally obtain 

In. „, L 



CeCs (<5o + l/6o) > Ao/2, 



An 



□ 



Reax{u,u) > -|1 A?/||^ + -|| Ar«||^,r + yll^lllr- 

By setting C5 = min {i, |, ^} and using ([21]l the proof is complete. 

Remark 1. It is clear from the proof that Aq > 4/1 and C5 < Z/4, so that Aq +0 
and C5 -> as Z — > 0+. This instability property will be confirmed in Remark [2j 

We can now give the 



Proof of Theorem [3l By Lemma [31 Problem © can be equivalently written as 
(|26p . The sesquilinear form ax in V is trivially continuous and, by Lemma 21 it is 
also coercive when ReA > Aq. We then apply Lax-Milgram Theorem (see [11] p. 
376]) to get the existence of a unique solution m of ([9]) in V. By Lemma El we also 
have u G H^. 

We now suppose that h g H", s > 1. To recognize that u E H'^^^ we apply the 
same bootstrap procedure applied in Lemma El More precisely, we shall prove that 
for any n G No we have 

(52) ue ^j^^ U|r e 77'"'"^"+2,n+4}(p)^ 

from which our claim follows for n large enough. We prove ()52p by induction on n. 
To prove that ([5^ holds when n = we recognize that, by ©i, 

so by elliptic regularity (see [26]) we have u e iy™in{«+2>7/2} j^j-j^^ which is the 
required regularity on f2 when n = 0. By the Trace Theorem we then have G 
^min{s+i/2,2}(-p)^ Hcncc by (111)2 wc have -Apu + M|r e i?™^'''^} (p). By the 
isomorphism property of — Ar + 1 wc then get U|r G ^min{s+2,4} which completes 
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the proof when n = 0. To complete the induction process we now suppose that 
([52| holds. Arguing as in the case n = by ((9])i we get 

Au = \u-he i7«"n{*'."+7/2} (f^) and U|r G iJ"iin{«+2:4+»} (D 

so by elliptic regularity u G ijmin{s+2,n+9/2}j-^-|^ gy ^^le Trace Theorem we then 
have € iJ™{«+i/2,n+3}(p)^ ^y using ([91)2, -Aru + U|r e H«"n{s:»+3} (p)^ 
before U|r € 7J™™{s+2,"+5}(p-)^ completing the induction process. 
Finally, to prove ([24|) we set up the operator Ax : £'(ylA) — > i?' , where 
D{Ax) ^{{u,v) e : (AM)|r = ku^ + lArv} 

and 

V^/ - IAtv + \v)- 

One easily sees that D{A\) is closed in H^^^ , so it is an Hilbcrt space with respect 
to the scalar product inherited by it. Moreover A\ is bounded, and u € H'^^^ solves 
^ if and only if m S D{Ax) and A\u = h. By previous analysis Ax is bijective, so 
(p4)) follows by the Closed Graph Theorem. □ 

4. Analysis of Problem ([T]) 

We will use here the results of the previous section to analyze Problem ([T|), thus 
proving Theorems [1] and [S] We start by setting up the unbounded operator A : 
D{A) CH Hhy 

(53) D{A) ={(u,w) e : (Au)|r = ku^ + lArv} 

and 

Our main results are a consequence of the following one. 

Theorem 4. Operator A generates an analytic semigroup {S(t),t > 0} in H, and 

(55) \\S{t)\\ciH)<e^°' forallt>Q, 

where Aq is the positive number given in Theorem so {S{t),t > 0} is quasi- 
contractive. 

Proof. We introduce the unbounded operator B in H by D{B) = D{A) and B = 
A — XqI. Then, given any u € D{B), we have that u solves ^ when A = Aq and 
h = -Bu. Hence, by ([26)l . 

(56) {Bu,u)h = —{h,u)H = —a,Xa{u,u) for all u e D{B). 

Then, by Lemma |4] we get that 'Re{Bu,u)H < 0, for all u E D{B), i.e. S is a 
dissipative operator in H. Moreover, by Theorem [3l R{I — B) = H . We then 
apply [28l Theorem 4.6, p. 16] to get that D{B) is dense in H. Moreover, given 
any u G D{B), hy Lemma |4] and ((56|) we have 

(57) Rc{-Bu,u)h = Rcaxo{u,u) > C5\\u\\y, 
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while by (|56|) and the continuity ol a\„ 

(58) |Im {-Bu,u)h\ < |aAo(u,^^)| < CMv 

for some Cg = Cg{k,l, N,n) > 0. Combining (|57p and ([58| we get that — i? is a 
densely defined m-sectorial operator in H. We then apply semigroup theory (see 
for example [21[ Theorem 5.9, p. 37]) whieh shows that B generates an analytic 
contraction semigroup {T{t), i > 0} in iJ, and consequently A generates an analytic 
semigroup {S{t), t > 0}, given by S{t) = e^o*T(i), i > 0, so clearly ^ follows. □ 



Now we can give the proofs of our main results. 



Proof of Theorem [H By Theorem [4] the operator A generates the analytic, and 
hence differcntiable, quasi-contractive semigroup {S{t),t > 0} in H. Then, by 
semigroup theory (see [28l §4.1]) given any uq E H there is a unique solution 

(59) ue C{[0,oo);H)nC\iO,oo);H) 



of the abstract Cauchy problem 
(60) 



u'{t) = Au{t), t > 

u(0) ^ Uq. 



Clearly, (|59p is nothing but (|52|) . and (|60|) is the abstract form of Problem ((T|). 
Moreover, ^ is nothing but ([55]) due to Lemma [31 Next, by using the differ- 
entiability property of the semigroup {S{t),t > 0} and [28l §2.4] we get that 
u e C°°{{0,oo);H) and consequently Bu = Au- Xqu = u' - Xqu e C°°{{0,oo); H). 
By ^ (when s = 1) then we get that u g C°°{{0,oo); H^). A standard boot- 
strap procedure then gives that u £ C°°{{0, oo); i/^"+^) for all n E N. By Morrey's 
theorem (see for example [51 Corollaire IX. 13] we then get that ([S]) holds. □ 

Proof of Theorem [2l We introduce, by recurrence on n € N, the space 
(61) L»(B") = {ue D{B''-^) : Bu G D{B''-^)) 

endowed with the graph norm 



n 



i=0 

By Theorem [3] it is immediate to recognize that 

(62) D(S") = {u e . (A*u)|j, ^ A:(A*-iu), + ;Ar(A'-iu)|r), t = 1, . . . , n}. 

and that the graph norm it is equivalent to the norm of i/^n+i jn^j-oduced in 
Section [3] Since i? is a dissipative operator in H and R{I — B) = H we are able 
to apply the procedure outlined in the proof of [5j Theoreme VII. 5] (see also [6l 
Chapter 1]) in the real case, which works as well in the complex one. Consequently, 
since uq e _D(i3"), we get 

u £ C([0, oo); D{B'')) D C\[Q, oo); £i(B"-i)) n . . . C"([0, oo); i?) 

which, by previous remark, is nothing but ([7]). Finally, if uq S C°°{^1) and ([6|) holds 
for all i G N we apply previous analysis, for any n £ N, together with Morrey's 
theorem to get (|8]). □ 



14 



JUAN LUIS VAZQUEZ AND ENZO VITILLARO 



5. Limit behavior as Z — >■ 0+ 

This section is devoted to study the limit behavior of the solution of Problem ([T|) 
when Z — > 0^ . The limit problem, at least formally, is given by 

{ut — Am = in Q = (0, oo) x il, 

ut ~ kuu on [0, oo) X F, 

u(0,x) = ui){x) in ri, 

which has been studied in [32] (see also [3]). We recall the following definition and 
result from that paper. In what follows we restrict to the real-valued case. 

Definition 1 ([32] Definition 1]). Given uq e H^{il) we say that 

(64) ueC{[0,T);H\n)), T > 0, 
is a weak solution of ([63]) if 

(65) u(0) = uo 
and 

T r 1 '■-^ 



Jn Jo Jn ^ JQ Jr 



(66) - / uipt+ I I VuVip + - / I uipt = 

Jo Jn 

for ah if e C^{{0,T) x IT). 



Theorem 5 First part of Theorem 1]). If N > 2 there is uq S C°°{n) 

satisfying the compatibility conditions 

A"uo = k{A"-^uo)^ on T for all neN, 

such that Problem (1631) has no weak solutions. 



The first step in our analysis is the following 
Lemma 5. Theorem\^ holds also if (|64|) is weakened to 
(67) u€C^{[Q,T)-H\^)), 
that is it concerns also weakly continuous solutions. 

Proof. Looking at the proof in the quoted paper one immediately sees that the 
continuity of u was used only at two places: at first in order that (|66p makes sense, 
and at second to recognize that the functions t (u(t),$^) are continuous in 
[0,T), where (•, •) denotes an equivalent scalar product in H^{0) and n e N 
are eigenfunctions of a suitable eigenvalue problem, which belong to C°°(ri). Both 
facts continue to hold when (l64l) is weakened to (l67l). □ 



We can now state the main result of this section. 

Theorem 6. Let Uq € H^(fl) be an initial datum such that Problem (j63p has no 
weak solutions u S Cu,([0, T); iJ^(r2)) for any T > 0, and denote by the solution 
of ^ corresponding to uq and I given by Theorem]^ Then, for any T > 0, we 
have 

niax^ ||u'(i)|| jLfijs^-) — > oo as /— >0+. 



^which is real valued since uq is real valued 
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Proof. We suppose by contradiction that (|68|) fails for some fixed T > 0. Then 
there is a sequence In — > 0+ such that, 

(69) ||u"||c([o,T];Hi(f2)) < Cio forallneN, 

were we denoted u" = u'" for simphcity, and Cio = Cio{T,uo,n) > 0. Since, by 
Theorem [U we have m" € C°°((0, oo) x Q) we are allowed, for any t e (0, T), to 
multiply the heat equation by a test function v € (17) and to integrate by parts 
to get 



(70) 



< |lu"(t)|Ui(,,)|lVH| 



By a standard density argument (|70|) holds true for all v G HQ{^l) and then, by 
we get the second estimate we need, that is 

(71) IKIU-((0,T);ff-i(O)) < Clo for all neN. 
By ((69|) we get that, up to a subsequence, 

(72) weakly* in L°°((0,r);i/i(f7)). 

Moreover, by combining (j69|) with (j7ip and using the compactness of the embedding 
H^{n) ^ L-^(f7) ^ iJ^^(ll) we also get, by Aubin-Lions compactness lemma, that 

(73) u''^v strongly in C([0,r];iJ"i(f7)). 

Moreover, by ^ we get that u" u weakly* in L°°((0, T); iJ^^ (!])), while by 
((75)) li" ^ u in the same sense, so u~ v. Then we can combine (|72P " ([75)) to 

(74) u" ^ u weakly* in L°°{{0,T); H\n)) and strongly in C{[0,T]; R-^n)). 

We now claim that u is a weak solution of in the class (|57)) . Once this claim 
is proved the proof is complete, since we are in contradiction with Theorem [5] as 
extended by Lemma [5l By applying [29l Theorem 2.1], we get 

(75) u€Cu,{[0,T];H'{n)). 

Moreover, by ([71)1 it immediately follows that (p5|) holds. Multiplying the heat 
equation in ([T]) by a test function ip £ C^{0,T) x il, integrating by parts in ft, 
using the boundary condition in ([1]), using (jlOp and finally integrating by parts in 
time twice we get that u" satisfies the distribution identity 



(76) - 



JO Jn Jo Jn Jo Jt Jo Jv 



for aU Lf € C;?°((0,T) x VL). By we can pass to the limit as n ^ oo in ([TS]). 
and we get that u satisfies (j66[) , completing the proof of our claim. □ 



Remark 2. Theorem [5] shows that the instability property of Ag = Ao(/) pointed out 
in Remark [T] does not depend on our estimates. Indeed, if one could improve them 
to get a sequence Z„ 0+ and Xo{ln) < A < oo, then by @ one would contradict 
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6. Open problems and final remarks 



Although Theorems [THl] give existence and uniqueness of solutions to Problem ([T]) 
in a Hilbert framework, building a satisfactory theory for C°° initial data, many 
interesting problems are still open, both of theoretical and applied nature. 

1. We are not able to produce a satisfactory regularity theory in even order 
spaces 7?^", n > 1, which is particularly bad for n = 1. A new estimate in 
V would be necessary. 

2. The extension of the analysis to more general problems, like the ones con- 
sidered by [9]-[10] has still to be done. In particular. Lemma [3] has to be 
properly extended. 

3. Our arguments, which are based on Lax-Milgram theorem, cannot be ex- 
tended to the case of Banach spaces. Now, it would be natural to consider 
the case of uq G W^^P^n), Mo|r e W^'P(T). 

4. We do not give explicit representation formulas of the solution u, even for 
regular data. In particular we are not able to apply the Fourier method, 
which would be based on the study of the eigenvalue problem 



The elliptic theory developed in Section [3] allows to prove in a simple way 
that 

S := {A e C : ((77)) has a nontrivial solution u &V} 

is at most countable, but it is far from giving an exhaustive spectral theory, 
since the operator A or equivalently Ax is not symmetric in H. Actually 
formula suggest suggest some symmetry of the operator A\, but in a 
framework of Krein spaces. This study is left to specialists in Krein spaces 
theory. 
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(77) 




in ri, 

on r. 
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